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Abstract 

The entropy production rate for an open quantum system with a classically chaotic limit has 
been previously argued to be independent of h and D, the parameter denoting coupling to the 
environment, and to be equal to the sum of generalized Lyapunov exponents, with these results 
applying in the near-classical regime. We present results for a specific system going well beyond 
earlier work, considering how these dynamics are altered for the Duffing problem by changing h, D 
and show that the entropy dynamics have a transition from classical to quantum behavior that 
scales, at least for a finite time, as a function of fi^/D. 
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Consider a quantum system with a nonlinear classical limit: Non-classical effects depend 
on the size of Planck's constant h compared to the characteristic action. Further, the 
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system-environment interaction as measured through some parameter D, is crucialpj. The 
dynamics of the classical limit of the problem are important [2] particularly through the 
classical Lyapunov exponents A. It has recently been found that the quantum entanglement 
rate for chaotic systems shows a valuable speed-up [3l but this is to be balanced against 
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the observed enhanced decoherence effects for chaotic systems in the classical limit |5|, 16| . 
However, there is increased stability against fidelity decay deep in the qiiantum parameter 
regime leading to the proposal to 'chaoticize' quantum computation [8]. This complex 
multi-parameter quantum-classical transition is fundamental, poorly understood, and also 
valuable in understanding the behavior of quantum devices. 

A recent analysislol], summarized below, suggested that headway could be made in char- 
acterizing the full range of behavior by considering composite parameters and scaling. That 
is, the quantum-classical difference as measured by some quantity QCd{h, D, A) should be 
the simpler function QC'^(C) of a single composite parameter C = h°'D^X'^. Evidence has 
begun to accumulate lOj] supporting this perspective. These come mostly from studying the 
effect of changing D, h on time-independent (usually from t oo) measures QCd- The 
change with A is harder to study since the classical phase-space changes along with A. A 
different but related issue is the non-equilibrium statistical mechanics of a nonlinear quan- 
tum system as measured through the system's entropy dynamics. A powerful result of broad 
interest is that the entropy production rate for an open quantum system with a classically 
chaotic limit is independent of h and D and is equal to the sum of generalized Lyapunov 
exponents However, this has been verified only in the classical limit, and despite 

the considerable interest in this, there are few useful results away from this limit. 

In this Letter, we start with the argument that quantum-classical distance can be mea- 
sured sensibly with a quantum system's linear entropy. We then study the entropy dynamics 
for the chaotic Duffing oscillator as a function of h, D to obtain several novel results that 
considerably extend results on entropy decay as well as generalize the scaling results. Specif- 
ically, the Lyapunov exponent dependence is shown to be valid only for a small parameter 
range and for times. We look, more usefully, at the time- dependent entropy itself which 
unexpectedly shows scaling with a singleparameter = fi^/D, thus generalizing previous 
results from time-independent measuresjol, [l^. That is, behavior from widely varied h,D 
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collapse onto curves that depend only on (q, which we explain on the basis of an expansion 
in (q, as well as direct comparison of dynamics. This enables the characterization of entropy 
dynamics over a much wider range of parameters and times than previously attempted. We 
show dynamical regimes which we term (I) classical, (II) semi-classical and (III) quantum, 
with a smooth transtion between these regimes with increasing (q. 

We begin with the Master equation for the evolution of a quantum Wigner quasi- 
probability pw under Hamiltonian flow with potential V{q) while coupled to an external 
environment 

dpw 



dt 



L, + Lg + T (1) 
= {H,pw} 

;,2»(_l)n g2n+l^(g) g^n+l^^ 

+ 2^dp{ppw) + Ddlpw (2) 

The first term, the Poisson bracket Lc, generates the classical evolution for pw- The terms 
in h, are the quantal 'correction' terms (denoted Lq). The environmental coupling (T) is 
modelled by the diffusive D term and the dissipative 7 term. We assume, as typical, short 
time-scales or high temperatures such that the 7 term is negligible. A QCd can then be 
considered by propagating the same initial condition with Lq + Lc + T, compared to using 
only Lc, or more appropriately using Lc + T from above. 

If QCd is the difference between the expectation values of an observable, it becomes 
strongly dependent on the observable. For example, even when the centroids of a quantum 
and classical distribution are behaving identically, differences exist in higher-order moments. 
Further, measures such as the time when the QCd hits a pre-defined value introduce sub- 
jectivity. Moreover, while powerful in the abstract, it is inherently unphysical to propagate 
something both classically and quantally. Some of these problems can be avoided by monitor- 
ing the quantum entropy, which does not measure distances but directly addresses relevant 
issues of information. The linear or Renyi entropy of second order 5*2 is also the natural 
logarithm of the purity P as 5*2 = ln(P) = ln[Tr{p^}]. Note that P = 2TThTr{p'^} where 
the Tr now represents integration over all phase-space variables. This has been extensively 
studied and for a system with a classically chaotic limit, it has been argued ^, ^ that in the 
weak-noise, small-h classical limit, —dS2/dt equals the sum of the positive classical Lyapunov 
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exponents.^ More careful considerations generalize this to a weighted sum over Lyapunov 

assical limit itself, this should arguably be further generalized to 



exponents [sl- For the c. 
time-dependent versions 



13!]. That is, although the previous results apply in some limits or 



special cases, even the classical behavior is not fully understood. Less is known about the 
quantum system, particularly the impact of changing scale or noise through h, D, which is 

what we address below. We work with the Hamiltonian H = ^ Bx^ + + Ax cosiut). 

This is the Duffing oscillator, which as a 1-dimensional driven problem with a quartic non- 
linearity is one of the simplest flows with a rich phase-space structure and hence is a paradig- 
matic problems in Hamiltonian chaos. The quantum version has also been frequently stud- 
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ied, including for decoherence issues [61,11^. We briefly review the behavior of the classical 
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FIG. 1: Entropy production rate 52 (where 52 = ln(P) for states with purity P) for the quantum 
Duffing oscillator with m = 1, B = 10, C = 1,^ = l,a; = 5.35 and h,D as indicated, showing 
a wide variation in behavior. The initial conditions are Gaussians in the chaotic region with 
{x)t=o = 1.0, the spread cr^ = 0.05, {p)t=o = 0.0 and the spread dp is set by the constraint of 
defining a minimum-uncertainty state and hence by the particular value of h. 

density pc in the limit of only the Lc + T evolution. As a result of chaos due to Lc alone, 
pc increases fine-scale structure exponentially rapidly, with a rate given by a generalized 
Lyapunov exponent. When the structure gets to sufficiently fine scales, the noise T becomes 
important, and it acts to decrease, or coarse-grain, fine-scale structure. These competing 
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effects can be profitably studied using the measure 

2 ^ Tr[peVVj ^ Tr[|Vp,p] 
^ - Tr[(p,)2] Tr[(p,)2] ^ ^ 

whence is approximately the mean-square radius of the Fourier transform of p, mea- 



suring the structure in the distribution [ij]. Most importantly, Eq. ([2]) yields the identity 
dS2/dt = -2Dx'fl Q 

with this valid classically or quantum mechanically, that is, with 
both 5*2, computed for pc or respectively. For a classically uniformly chaotic system, 
the dynamics of can be written approximately 15|] as a competition between chaos and 
diffusion as 

^ ^ 2Ax' - ADx'. (4) 

This implies that settles after a transient to the metastable (that is, constant for finite- 
time) value 

X'* = A/2D (5) 



where A is a p dependent generalized Lyapunov exponent [15|, [ITJ • This classical argument 
leads to the argument 4, ^, ^ that quantum entropy-production rates are equal to gener- 
alized Lyapunov exponents. This applies to a greater range of parameters than might be 
anticipated because decoherence suppresses quantum effects. While this behavior has been 
shown in several instances, it does not capture the complete picture, particularly the effect 
of changing h,D. We show this in Fig. ([T]) plotting dS2/dt for the Duffing problem with 
m = 1, B = 10,C = 1, A = l,uj = 5.35, as previously usedGj. The behavior, over a wide 
parameter and time range is quite complicated. If a subset (all of those with h = 0.1) 
are plotted for a short time (t < 15) as in[6], they show the classical Lyapunov exponent 
entropy-production behavior . This is valid only for some small range of parameters 

and short times. There has been a suggestion of a superposition of classical and quantal 
exponential decay HI] for the purity. This would lead to a crossover transition within a fairly 
narrow range from one constant value to another in Fig. (1), which we do not see. Other 
ways of considering the data (as in Fig. (2) below) also do not support this. In general the 
search for these small regimes of linear decay for entropy is not as helpful as understanding 
the broader parameter dependence. 

To do this, consider as in Fig. (2), Tr{p^(t)}/Tr{p^(0)}. Since the y axis is logarithmic, 
we are effectively looking at ln(Tr{p^(t)}) - ln(Tr{p^(0)}) = S2(t) - ^2(0) = S{t) for our 
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pure-state Gaussians. This shows useful organization invisible in Fig. (1) due to the small- 
scale variation in a narrow range. Most interestingly, the entropy dynamics for the wide 
variety of parameters considered is captured entirely for the times shown by the composite 
parameter fi^ / D = Co, even though a wide range of behavior, not obviously characterized as 
exponential decay, is seen as Co is varied. Larger Co corresponds to high h or low noise D 
or both, and remains closer to a pure quantum state for longer times, which makes physical 
sense. Note also that there is some Co dependence for the time-scale of scaling, with a 
long-term separation of curves. 
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FIG. 2: Evolution of the normalized purity for the same states as Fig. (1) in the quantum Duffing 
oscillator. Scaling is observed relative to the parameter Co = h^/D. (I)'Classical': Co = 2 a)?i = 
0.1, D = 5 X 10"^ h)h = 0.2, D = 2 x 10-2, c)h = 0.5, D = 0.125, (II) 'Semi-classical': Co = 40 
d)h = 0.004^2, D = 10-4, e)h = 0.1, D = 2.5 x 10^^, l)h = 0.2, D = 10~^, (III)'Quantum': 
Co = 100 g)h = 0.1, D = 10-'^, h)h = 0.5, L> = 2.5 x lO^^ i)^ = 1,D = lO^^. 

We understand this Co dependence by considering quantum corrections to the classical 
dynamics, which depend (see Eq. (|5])) on the derivatives of the Wigner function. Given that 
the second derivatives d'^pw oc x^, these corrections scale as 

r _ fc2n ^^"^^'^(g) ^^"^V^ _ ^2n 2n+lyi2n+l)( ^ 

^9 ~ Qq2n+1 Qp2n+1 ~ ^ X V [X), (DJ 

where V^^^ denotes the rth derivative of V. When the phase-space distribution hits a 
metastable state such that settles to the fixed value A/2D, the difference between the 



quantum and classical evolution may be estimated to depend on 

^ = ;^2n^"+l/2£)-(n+l/2)y(2n+l)^^^^| 



where, since x is a 'length' in Fourier space, we have that a; ~ x y'^^/A. This is 
essentially the same 

and is also the root of the suggestion in Ref. to search for scaling. Therefore, the first order 
quantum corrections in a semi-classical regime should scale, in complete generality, with the 
single parameter (. The particular form of ( is decided by the details of the Hamiltonian 
and the difference between the quantal and classical propagators. For the Duffing problem, 
the only quantum term of Eq. ([6]) comes from the 3rd derivative of the quartic term whence 
Eq- © gives that the quantum term goes as C = ^^X^j fo^' other form of the potential, 
we expect different corrections and hence different scaling as below. 

We now use this in an expansion technique for entropy dynamics that may be applied 
in general. In the Duffing problem, even though 5*2 is not a simple function of (, a scaling 
relationship still obtains in the two parameters h, D as follows. To zeroth order, the classical 
and quantal phase-space distributions are the same, pwo ~ Pc, and x^o = Xc; where the 
entropy production rate S2qo = ~'^DXqo and the numerical subscripts on Xq, Pw, S2q indicates 
the order of the approximation. We now use the results from Eq. ( pfTI) that the quantum- 
classical distance for this system behaves as h'^x'^- To first order we insert the zeroth order 
solution in this to write 

Pwi ~ Pwo + afi^X^Pwo = Pc + ah^x^Pc (8) 

where a is constant for the meta-stable state, but time-dependent in general. We substitute 
this in Eq. ([3]) to get that Xqi ~ Xc + ^^^xt- Corrections from the denominator of Eq. ([3]) 
are of higher order, and also tend to cancel the higher order corrections from the numerator. 
We insert this first order quantally corrected form for the dynamical term into Eq. (jlj) to 
get that to first order in fi^^ x^ obeys 
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dt 

and in parallel to Eq. ([5]) we get that 



2A(x^ + ah\') - 4Dx' (9) 



leading finally to 52^1 = —2DXgi = — A(l + ^) that is, the quantum correction scales 
as fi^/D. This expansion around the metastable state can occur only when the growth of 
structure is balanced by noise, only when is large enough that the diffusion term becomes 
relevant. Since a is in general time-dependent, at each value of Co we expect a different 
entropy dynamics, as in fact we see. In sum, this expansion for the entropy dynamics around 
the metastable state yields a h^/D dependence for entropy, although the time-dependence 
itself is not easy to extract. 

This expansion must fail for arbitrarily large (o, in the quantum regime. Here an alternate 
approach applies: the Poisson bracket term is neglected and the dynamics are given approx- 
imately by the competition between the Lg and T terms alone. To compare them, consider 
Lg-. For the Duffing system, there is a third-derivative of pw multiplied by x (resulting in 
this acting like a 2nd derivative overall), compared to the 2nd derivative from the diffusion 
term lid]. This means that the terms have essentially the same scale, with quantum dynam- 
ics continuing to add structure and the noise smoothing it out. The entropy-production 
then depends only on the ratio of the parameters multiplying these terms which is again 
h'^/D = (o. This last parameter regime is consistent with recent resultsj?!, et al. Fi- 
nally, consider some details of the time-dependence: The rate of purity decay decreases with 
Co- Physically, the time- asymptotic dynamics are dominated by essentially classical diffu- 
sive behavior, with a common final state (the natural invariant measure) for all p^. Since 
Tr{p^(0)} oc h^^ (see above), the time-asymptotic value of Tr{p^(i(:)}/Tr{p^(0)} oc h~^. 
With the different rates of purity decay, the system approaches the time-asymptotic state 
later as Co increases. Further, within each Co, the different values of h separate out from the 
scaling curve as the final diffusive regime kicks in, as seen in Fig. (2). 

The values of Co where these regimes change is in general determined by the parameters 
of the potential, i.e. by the quantity labeled as a in Eq. ([8]) above. Given the continuous 
behavior as a function of Co the actual transition is subjective. In Fig. (2) we label what cor- 
responds to rapidly decohering and hence essentially classical behavior as (I), the relatively 
slowly decohering and hence deep quantum behavior as (III) and in-between 'semi-classical' 
behavior as (II) in the three sets of curves with (Co = 2, 100, 40) respectively. That is, for this 
potential, empirically Co = Cc ~ 10 sets the approximate upper limit of the rapidly decoher- 
ing regime (I), and by extension the quantum regime (III) kicks in at Co = Cg ~ Cc ~ 100- 
note the same scaling also holds (results not shown) for other diagnostics as well as for very 
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different parameters for the Duffing oscillator, A = 10,0; = 6.07, a regime of significantly 
increased chaos j^. 

In conclusion, our results strengthen the argument that it is valuable to study the behavior 
of nonlinear open quantum systems through the scaling behavior of appropriate diagnostics, 
as recently suggested^]. In particular, this is used to study the non-equilibrium statistical 
mechanics of an open quantum system with a classically chaotic counterpart over a wide 
parameter range in h, D. We show that the entropy dynamics of this system can be dra- 
matically different from the broadly-accepted Lyapunov exponent dependence which is only 
valid in the classical limit (and is itself arguably suspect We show a h^/D scaling in the 
time-dependent entropy dynamics, although the particular form of the scaling is expected 
to depend on the form of the nonlinearity in general. 
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